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^ ; Abstract 

A complete discussion of the constraints on the Michel parameters and the am- 
biguities of their interpretation is presented. Estimators of new physics, optimized 
for a very wide class of hypotheses and models are proposed. 



a,: 
a: 

1 Introduction 

X 

The leptonic decay of a charged lepton l x — > l 2 U\ i>2 is described by its Michel parame- 
ters [1-4] Complementary parameters must be introduced to describe the secondary 
lepton polarization. The measurements of all the parameters and the cross-section for 
the V\ I2 — > l\ V2 reaction allow to put upper limits on all the non-V-A amplitudes M, M 
and "prove experimentally" the V-A structure of the weak interaction. 

In the case of the \i lepton, this program was followed and successfully completed 
fifteen years ago || [7|. 

For the r lepton, it is customary [7-9] to present the experimental results, in a sim- 
ilar way, as upper values of the possible coupling constants. However, the validity of 
the Standard Model is today so well established, that the interest of r Michel parameter 
measurements is no more to prove the V-A structure but to look for small deviations from 
its predictions that would be evidence of new physics. Moreover, the additional measure- 
ments, which are necessary to complete the experimental proof of the V-A structure, will 
not be performed in a foreseeable future, while large statistics of r pairs will soon be 
available at B-factories, allowing more precise measurements of the Michel parameters. 
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Therefore it seems worthwhile to look for the combinations of the measured parame- 
ters which are the most sensitive to non-standard effects. This is the aim of the present 
paper. 



2 The general framework 

If only the momentum of the final state charged lepton is measured, the decay of a 
polarized r is entirely described, in its centre-of-mass, by the distribution 

where P T is the r polarization, 6 the angle between the polarization and the charged 
lepton momentum and x = E\j '_E , ™ ax the normalized lepton energy. 

Taking advantage of the weak interaction short range, the decay can be represented 
by the most general four fermion contact interaction [1-4], written below in the helicity 
projection formalism [6,10-12]: 

- M = 4 % E alAUn W<(*r)jr 7 |r,), (2) 

j=S,V,T , 

i,j=R,L I = e,fJL . 

G T i is the absolute coupling strength; the gj, are ten complex coupling constants de- 
scribing the relative contribution of scalar {V s = 1), vector (T = 7 M ), and tensor 
(r T = -^c^ u ) interactions, respectively, for given chiralities j, i of the r and the charged 
decay lepton. The neutrino chiralities n and m are uniquely defined for a given set 
{7,i,j} by the chirality selections rules: conservation for vector coupling, reversal for 
scalar and tensor. 



Table 1: Contributions of neutral currents to the gj~ coupling constants 



Neutral 
current 






Charged coupling 


constants 




V 


9ll 


9 V RR 


9 S LL 


9rr 






S 


9 V LR 


9 V RL 


9lr 


T 

9lr 


9rl 


T 

9rl 






9 S LR = 


ig T L R 


9rl = 


ig T R L 


T 






9lr 


T 

9lr 


9rl 


T 

9rl 








9lr = 


Sg T LR 


9rl = 


-§9rl 



The same matrix element describes also the contributions of possible lepton-number- 

The relationship between the neutral currents and the g] j 



violating neutral currents [[13 



constants, as given by the Fierz transformation, is summarized in Table |I| 

The matrix element (fj) can be used to compute the functions Wq and W\ of Eq. [[]. 
They are described by the four real Michel parameters: p and r\ for Wq, £ and £5 for W\. 




Restricted domain 



Figure 1: The allowed domain in the space of the three parameters, p, £, and ^5 

A convenient way to express the relation between the Michel parameters and the gjj 
coupling constants is to introduce jTJ| the six positive parameters 



(3) 
(4) 
(5) 





= \9 V R L \ 2 + \9 S RL + §9 T RL 


2 /16 


a~ 


= \9 v lr\ 2 + \9 s lr + 59 t l R 


2 /16 


P + 


= \9U 2 + \9rr\ 2 /4, 




P~ 


= \9 V LL \ 2 + \9 S LL \ 2 /^ 




7 + 


= (3/16)\g s RL -2gl L \ 2 




7~ 


= (3/16)\g s LR -2gl R \ 2 





They satisfy the relation 



1 = a + + cT + (3 + + /T +7+ + 7" (6) 
= \9 V rr\ 2 + \9 V L r\ 2 + \9 V r L \ 2 +\9 V LL \ 2 

+\ (\9 S rr\ 2 + \9 S LR \ 2 + \9k\ 2 + \9ls\ 2 ) + 3 (\9Ir\ 2 + \ 9 T r L \ 2 ) , 
which means that the normalization is absorbed in the definition of G r \. 
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The p, £, and £5 parameters are given by 

p = \^ + + n + {i + + r) , 

£ = 3(«--a + ) + (/r-/3 + ) + ^(7 + -7") , (7) 

^ = + -r) ■ 

In geometrical terms, the point of coordinates p, £, and £5, in the space of the pa- 
rameters, is the barycentre of six points, A , B , and C ± with the weights a , f3 , and 
7 respectively. Since the point B~ lies on the A + C + segment, and B + on A~C~, the 
allowed domain is just the tetrahedron A + A~C + C~ (Fig. [I]). 



Table 2: Values of the coupling constants in the Standard Model and upper values compatible 



model prediction: £ = 


= 1, 




= 3/4 
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p = £5 = 3/4, £ = 


1 












1 T l 

\9lr\ 


\9rl\ 
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1/V3 


1/V3 




SM 














maximum 

p = £5 = 3/4, f = 


1 







1/2 





The Standard Model prediction, |g^ L | = 1, implies p = £<5 = 3/4,£ = l and is 
represented geometrically by the point B~ , but /3~ = 1 does not imply \g][ L \ = 1, and 
the location of B~ on the A + C + segment introduces further ambiguities. 

The value (3~ = 1 is also obtained if = 2 and all the other constants equal 

to zero, and there is another possibility to reproduce the parameters predicted by the 
Standard Model, namely (3~ = (3 + = a~ = -y~ = 0, 7 + = 3/4, and a + = 1/4. So, with 
£ = 1, p — £5 — 3/4, \gft L \ can reach an upper value of 1/2 and it is also possible that all 
the constants vanish but g^ L and g\ L . In this last case, we get 

|^ L | 2 + 12|^ L | 2 = 4, 2|^ L | = -|^ L |cos0 5T , (8) 
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where 4>st is the relative phase of the two amplitudes. Accordingly, the upper possible 
values of \g RL \ and \g RL \ are 2 and 1/2 respectively. 

The upper values of the constants, compatible with the Standard Model prediction, 
are given in Table 0. More stringent experimental limits |15[] are not the consequences of 
the data but of additional hypotheses or constraints. As the physically interesting region 
is the neighbourhood of \g\j\ = 1, no relevant bound on the non-standard r left-handed 
couplings can be extracted from the measurement of p, £ and £<5 without additional hy- 
potheses. 



Since the Standard Model prediction is represented by B~ which is located on A + C + , 
a convenient set of parameters is given by two equations of the A + C + line and a third 
variable which specify the position on the segment. Using the equations of the faces 
A + C + A~ and A + C + C~ gives expressions which are positive for points inside the tetra- 
hedron. They are more usefully combined into 

n = ^[l + f-y^]=/3 + + a-+7- (9) 
= \9 V rr\ 2 + \gU 2 + \ \9 S rr\ 2 + \ \9 S lr? + SlsLl 2 , 

s R = l[p-m = P + + l^ (io) 

= \9rr\ 2 + j I 9rr 1 2 + ^ \9lr ~ ^9lr\ 2 > 

which, owing to their normalization, can be interpreted as fractional contributions of r 
right-handed couplings to its leptonic partial width and used to bound the gJ R amplitudes. 
For the third parameter, £ or p can be chosen. 




Figure 2: The allowed domain in the a, (3 plane for a given value of p 



The last Michel parameter r\ can be written 

V = 2Re[g R \l(g s LR + 6g'E R )/4} + 2Re[g^ R (g s RL + 6gl L )/A} 

+ Re[gllg S RR /2] + Re[g v RR g s LL /2] . (11) 
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This expression implies the inequality 

M<a + ^<l-p + ~/Wp), (12) 

where, a = a + + a~ , [3 = [3 + + /3~, and /3 ma x(p) is the greatest value of /3, compatible 
with the value of p. Fig. || shows that /3 max is reached for 7 = 1 — a — (3 = when p < 3/4 
and for a = for p > 3/4, leading to the rather weak bound: 

M<l-|o (P<^)> N<2(l-p) (p>~) (13) 

3 The restricted domain 

If only charged vector currents are present (gfj = 0,gfj = 0), we have 7+ = 7" = and 
the allowed domain is the A + A~B + B~ tetrahedron that we call the restricted domain. 

It is interesting to remark that, for a very wide class of models and hypotheses, the 
allowed domain for the Michel parameters is this restricted domainQ and that getting a 
point outside it requires the conspiracy of scalar and tensor couplings fI3fl . 

We give below a list of general hypotheses leading to the restricted domain. We use 
for that the Lorentz covariance of the charged and neutral currents and the properties of 
the Fierz transformation only. More specific constraints associated with the hypothesis 
of a single boson exchange can be found in [12] and |L6 1 . 



When the allowed domain of the parameters is the restricted one, the point B~ , which 
is the Standard Model prediction, is one of its vertices, therefore there is a third positive 
quantity whose vanishing is associated with this point. We use for it 1 — 4p/3. Its precise 
physical meaning and the ambiguities in the interpretation of the measurements depend 
nevertheless on the hypothesis that lead to the restriction. 

9 s =0 

If all the scalar coupling constants vanish, the condition a + = implies 7 + = 0, 
therefore the complete domain is not allowed. Defining the new parametrization 

"ns = \9rl? , 7ns = 3|#L| 2 > "ns = l#L| 2 > 7ns = %L| 2 > ( 14 ) 
we get the relations 

1 = c4 s + c% s + I3 + + I3~ + 7ns + 7ns > ( 15 ) 

P = ^(/? + + /n + ^(7Ns + 7Ns) > 

£ = 3 (a NS - o4) + 09-- (3+) - jj (7+ s - 7 ns) , (16) 

£5 = 1(P--P + ) + \(^s-1ns)- 



1 This result is almost obvious from the "charge retention" formalism, where 7 + and 7 are related 
to (neutral) tensor currents Q] , and Table |l|. 
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The point of coordinates p, £, and £5 is now the barycentre of A , B , and two new 
points C^s = D ± whose coordinates are (1/4, =F5/3, ±1/4) but, since D + is located 
on the A + B~ segment and D~ on A~B + (Fig. [1]), the allowed region is the restricted 
domain. The third constraint can be written, in terms of the g^'s, as 

i-^ = bU + bLi 2 +2(bLl 2 + bLi 2 )- 

Measurements represented by B~ imply a V-A structure. The combination of charged, 
vector and tensor currents belongs to this class. 

If gfj — 0, we define 

+ _ I V |2 + 1 | 5 |2 - __ I V |2 - 1 | 5 |2 / 18 \ 

^nt _ \9rl\ j 7nt _ ^ > a NT _ Ifl'iiil > 7nt _ ^ \9lr\ j 

and obtain 

1 = a^ T + «nt + (3 + + f3~+ 7nt + 7nt , ( 19 ) 
3 

P = -^(P + + P~ + 7nt + 7nt) > 

e = 3 («nt - «nt) + (P~ - P + ) + (tSt - 7nt) > (20) 
= ^ (r - /? + + 7^t - 7nt) • 

Here the points C^ T and C^ T coincide with the points B~ and B + respectively. The 
allowed domain is then the restricted one. The third constraint is 

^-\p=\9 V RL \ 2 + \9 V LR \ 2 ■ (21) 

The vector couplings only are bounded by the measurement of p, due to the ambiguity 
between B ± and C^ T . 

The combination of charged, vector and scalar currents and neutral, vector currents 
belongs to this class. 

9rl = ^9rli 9lr = ^9lr 

The restriction of the domain is evident from Eq. [|. The third constraint can be 
written 

1 - \p = \9 V RL \ 2 + \9ln\ 2 + \ (\9l R \ 2 + \9 S RL \ 2 ) + 3 (lsLl 2 + I^LH • (22) 

Measurements represented by B~ imply a V-A structure. 

The combination of charged, vector currents and neutral, vector and scalar currents 
belongs to this class. 
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"V-A plus something"^ 



In an especially interesting family of models [ICj the decay is described by the ad- 
dition of a single, non-standard contribution to the Standard Model amplitude. Their 
predictions are perfectly transparent in the present geometrical presentation. They follow 
at once from the definitions of the parameters and the properties displayed in Table |l|. 

-If the non-standard contribution is a charged vector current, all the restricted domain 
is allowed. The only new prediction is rj = 0. 

-The hypothesis of an additional charged scalar current belongs to the second of the 
above defined classes with the further conditions a^ T = 0. The allowed domain is then 
the B + B~ segment (p = 5 = 3/4), since the points C^ T and B T are identical. 

-The contribution of a neutral vector is included in the same class with a^ T = 7^ T = 0. 
The allowed domain is again the B + B~ segment. 

-The hypothesis of an additional neutral scalar leads to 7 + = 7" = (third class) 
and (3 + = 0. The allowed domain is the two-dimensional one, spanned by the points 
A + , A~ , and B~ . The corresponding condition is p = £5 and the bound on r/ is stricter: 
\v\ < l-4p/3. 



4 Looking for new physics 

In the standard approach described in section [[J the indicators of new physics constructed 
with the Michel parameters are, besides the rj parameter itself, the two positive quantities 
V T R and S R , defined by Eqs. || and which bound the coupling constants and can be 
interpreted as (non-independent) contributions of new physics to the r decay. 

Using the world-average values of the parameters @, under the hypothesis of e-/i 
universality, yields 

V T R = 0.006 ± 0.028 , S T R = 0.001 ± 0.016 . 

If hypotheses are made that reduce the dimensionality of the domain, the determina- 
tion of the parameters can be improved by a constrained fit. For instance, if the domain 
is the B + B~ segment, there is only one free parameter. Neglecting the correlations be- 
tween the measurements and taking advantage of the near equality of the errors on £ and 
4£5/3, this parameter is merely their average £ and the quantity V R reduces to 

|[l-f] = |[l-|-|^] ■ (23) 

It is noteworthy that the same strategy can be followed under the much weaker 
hypotheses that imply the restricted domain. This is due to the fact that the measured 
value of p forces the point which represents the measurements in the parameter space to 
lie on the edge of the domain. 

Quantitatively, the difference of £ and 4££/3 is bounded by the inequality 

le - < 3(1 - |p) - (24) 



2 The title of this paragraph is borrowed from C. Nelson Ml 
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From the values of the errors |7J cr(£ — 4£<5/3) = 0.044 and 4cx(p) = 0.036, it is clear that 
there is no physically relevant information in £ — 4£<5/3 and that even the r-l universality 
prediction |12[ 5 = 3/4 is better tested by the measurement of p than by comparing £ 
and £5. 

Assuming only the restricted domain, the indicators of new physics are 

4 1 r z n 

1 - - p = 0.004 ±0.012 , -[1 -£] = 0.002 ±0.011 . 

o 

They bring a clear improvement of the sensitivity with respect to V T R and S T R . 

It must be noted that 1 — £ is not strictly positive but that excursions of £ beyond 1 
are severely limited by the inequality 

£-l<i(l-|p), (25) 
since cx(£) = 0.022 and 2<r(p)/3 = 0.006. 

5 Conclusion 

A complete study of the constraints on the Michel parameters and the ambiguities of 
their interpretation has been presented. 

It has been shown that, for a very wide class of hypotheses and models, which cause 
the same restriction of the parameter domain, the best indicators of new physics are the 
combinations, 1 - 4p/3 and (1 - 0/2 = (1/2 - £/4 - £5/3). 

Compared to the customary estimators, their sensitivities are roughly twice better. 
The third available parameter, £ — 4£<5/3, is better determined by the geometry of the 
domain and the value of p than by its measurement. 

A Appendix: Measurement of the parameters 

In the numerical exercise above, the error correlations were neglected and the e-p uni- 
versality was assumed. We will discuss briefly this two approximations. 

At low energy, where the r's are unpolarized, the p and r] parameters are determined 
by the single- lepton laboratory energy distributions and the £ and £5 parameters by 
the spin-correlated r + r~ decay distributions. Estimates of the covariance matrices for 
measurements at 4 GeV and 10 GeV can be found in ||. 

At the Z peak, the r polarization makes ambiguous the interpretation of a single tau 
leptonic-decay distribution and, since the transverse spin correlations depend on the Z 
couplings, only the helicity correlation of r + and r~, which is equal to -1, is used. If the 
decay distribution of a r~ in the channel a, is written, 

W a (x~) = f a (x-) + P T g a (x-) , (26) 

the correlated distribution for the a and b channels reads then 

W ab (x-,x + ) = f a {x-)f b {x + )+g a (x-)g b (x+) + P T [ f a (x~)g b (x + ) + f b (x + )g a ( X -) ] , (27) 
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where P T is the r~ polarization. 

For a hadronic decay, with the notations of [[HJ, x ± is the optimal variable u and the 
decay distribution is 

W(u) = f(u)[l+£ h P T u] , (28) 

where the ^ parameter is equal to 1 in the standard modelQ. 

For a leptonic decay, x ± is the normalized energy of the charged lepton, y = ( Ei/ Ef 13 ^) lab- 
Defining the parameters, 

P=1-~P, ~6 = Z-U6, V = —V, (29) 

3 3 m T 

which vanish in the Standard model, and the functions 

h (y) = 1(5 - 9y 2 + Ay 3 ) , h^y) = 1(1 - % 2 + 8y 3 ) , (30) 

Mv) = 1(1 - 12y + 27y 2 - 16y 3 ) , /i 3 (y) = 12(1 - y) 2 , (31) 

the decay distribution reads 

W(y) = f P M + p r9uM = T ^{ho(y) + [p + P^}h 1 (y)+SP T h 2 (y) + vh(y)} ■ (32) 

The presence of the polarization allows the measurement of a new parameter, 5P T , but, 
as previously mentioned, it also introduces an ambiguity in the interpretation of the first 
one which is now p + P T £. 

The same kind of ambiguity also arises in the e-p correlation. For P T = 0, only the 
parameters, 5^, 8 e , and the product are measurable ||. For P T ^ 0, neglecting fj 
for the sake of clarity, and keeping only terms of the first order in the Standard Model 
violating parameters, p, 5 and 1 — £, the correlated distribution can be written 

W(y e ,yn) ~ h (y^)h (y e ) - h^y^hifjje) + [pe + P<r€e]MVi*)MVe) ( 33 ) 

+[p M + PT€n]Mye)Myn) + + & + p ^pe + h)\My»)Mye) 

+S e h 2 (y e )[h 1 (y IJ ) + PtMVh)] + hMvJiMVe) + PrMVe)] ■ 

Even if P T is known, there is only three measurements to determine the four parameters 
Phi Pe, 6/u and £ e . The ambiguity is displaced but not suppressed^. 

Using several known values of the polarization, all the parameters can be determined 
from single-decay distributions. At the Z peak, and/or with polarized beam, the r 
polarization is a function of the production angle 9, hence the Michel parameters can be 
measured by the 9-y correlation [plj| . 

The other measurements [15,22-24] use all the hadron-hadron, lepton- lepton and 
hadron-lepton final states to obtain the complete set of parameters up to a global sign 
ambiguity which is solved, for instance, by the result of [21] . 

3 Since P T + = —P T - and = — £1 , the decay distribution is independent of the r charge. In 
the case of r — > is3ir, this last property is true only if no pseudoscalar variable constructed from the 7r 
momenta is used in the definition of u>. A more general analysis is presented in JlS] , 

4 Therefore, the analysis p0[ necessarily uses additional hypotheses. 
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To calculate the covariance matrix V of the Michel parameter measurements, we 
assume that, in an ideal experiment at the Z peak, similar to [15,22-24], all the decays 
into p, e, tt, p and a\ and their correlations are used to determine the Michel parameters, 
the hadronic, £/, parameters and the r polarization. Asymptotically V is given by 

where s labels one of the twenty classes of events, e-p, e-e, e-7r, etc, and N s is the number 
of events in the class. The computation is made straightforward by the fact that the 
distributions are quadratic functions of all the parameters, except rj which appears in the 
normalizations. 

Table 3: Ideal statistical errors on the Michel parameters (in %) and their correlation coeffi- 
cients for 2xl0 5 r + r~ pairs at the Z peak 





°{p) 




^ - 




2.3 


6.0 4.0 


7.3 


e 


1.3 


5.5 3.9 


7.2 






c(M) 


C{p,£S) 






0.24 


0.15 


0.18 


e 


-0.21 


-0.06 


0.10 



The computed covariance matrix is similar to its estimations || for measurements at 
lower energy. The largest correlation coefficients are C(p M ,^ /Lt ) = 0.82 and C(^,r]^) = 
0.42. The numerical values relevant for the analysis of Section ^ are given in Table [| 

For the r — > evv channel, the inequality o{$) > 3cr(p) is satisfied and the weights of £ 
and 4£5/3 in their optimal combination are 0.47 and 0.53 respectively. All the hypotheses 
made in Section || are verified. 

For the r — >• pvv channel, the weights are 0.43 and 0.57 but, owing to the correlations 
with rjfj,, the error on 5 is slightly smaller than 3<r(p). However, in a more realistic 
estimation, the inefficiencies in the identification of the various decay channels reduce 
the statistics for the classes of events with two analysed decays and increase it for the 
events with only one identified decay which contribute mainly to the measurement of p. 
Therefore, the same analysis scheme remains basically valid. 

From an experimental point of view, the universality hypothesis allows to constrain 



the value of tj^ by the measurement of the r — > evv and r — > pvv branching ratios [15, 20 



If the variation of the parameters is limited to the above defined restricted domain, 
the deviations of the parameters from their Standard Model values must have the same 
sign in the e and /i channels. Therefore the universality hypothesis can perhaps reduce 
the sensitivity to these deviations but complete cancelations are not possible. 
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